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Abstract. The field theoretic renormalization group is applied to the stochastic 
Navier-Stokes equation with the stirring force correlator of the form fc4-d-2e ^j^g j_ 
dimensional space, in connection with the problem of construction of the \/d expansion 
for the fully developed fluid turbulence beyond the scope of the standard e expansion. 
It is shown that in the large-d limit the number of the Feynman diagrams for the Green 
function (linear response function) decreases drastically, and the technique of their 
analytical calculation is developed. The main ingredients of the renormalization group 
approach - the renormalization constant, [3 function and the ultraviolet correction 
exponent w, are calculated to order (three-loop approximation). The two-point 
velocity- velocity correlation function, the Kolmogorov constant Ck in the spectrum of 
turbulent energy and the inertial-range skewness factor S are calculated in the large-d 
limit to third order of the e expansion. Surprisingly enough, our results for Ck are in 
a reasonable agreement with the existing experimental estimates. 
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1. Introduction 

One of the most interesting open problems inherited by the modern theoretical physics 
from the twentieth century is that of description of fully developed hydrodynamic 
turbulence on the basis of a microscopic model and within a consistent perturbation 
scheme [1]. A challenge which is still waiting to be answered is the derivation of 
anomalous scaling behaviour of the velocity correlation functions from first principles 
and calculation of the corresponding anomalous exponents within a regular perturbation 
theory, analogous to the e or expansions of the critical exponents in the theory of 
second-order phase transitions. 

The first difficulty is that the ordinary perturbation theory for the stirred 
(stochastic) Navier-Stokes (NS) equation (that is, the expansion in the nonlinearity) 
is in fact the expansion in the Reynolds number, a parameter which tends to infinity 
for the fully developed turbulence. Hence the necessity to rearrange (to sum up) the 
naive perturbation series. A similar problem is well known in the theory of critical 
behaviour, where it had been solved a long ago by means of the renormalization group 
(RG) and the self-consistency ("bootstrap") diagrammatic equations [21 [3]. The first 
way naturally leads to the famous e expansion (where e = 4 — d is the deviation of 
the spatial dimension d from its upper critical value d = A), the second leads to the 
alternative 1/N expansion (where N is the number of components of the corresponding 
order parameter). So far, however, those methods have had relatively limited success 
when applied to the problem of anomalous scaling in fiuid turbulence. 

An important difference is that the turbulence or, better to say, the stochastic NS 
equation, has no upper critical dimension, and the parameter e in the standard RG 
approach has completely different meaning. Namely, the correlation function of the 
random stirring force that provides the energy supply to the system is taken in the 
power-law form |1] 

(//) cx k'-"-'^, (1.1) 

where k is the wave number; more precisely, see [3], O [6] and section [2j 

The physical value e = 2 corresponds to the energy pumping by the largest scales, 
(//) oc 5(k), while d remains a free parameter and can be varied independently of e. 
The common point with the models of critical behaviour is that, in the both cases, the 
limit e — > corresponds to a logarithmic (exactly renormalizable) theory and e serves 
as the formal small expansion parameter in the RG approach. For this reason, we use 
the same symbol e for the stochastic NS problem (where it is also sometimes denoted 
as e = y/2). 

The results of the RG analysis of this model are reliable and internally consistent 
for asymptotically small e, while the possibility of their extrapolation to the physical 
value e = 2 and thus their relevance for the real fluid turbulence is far from obvious. Of 
course, the physical value of e = 4 — c? in the RG theory of critical phenomena is not 
small, either. But there, no qualitative changeover in the behaviour of the system is 
expected when e increases from the region £ <C 1 to real values e ~ 1, and the possibility 
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of this extrapolation is usually not disputed. The situation in the RG approach to the 
NS model (11 .11) appears different. New physical effects are encountered as e grows, and 
they can be easily lost or misrepresented if the e expansion is too naively applied. One of 
them is related to the so-called sweeping effects, the transport of small turbulent eddies 
as a whole by large-scale ones. The sweeping leads to strong dependence of the velocity 
correlation functions of the integral (external) turbulence scale L. The other effect is 
the crossover from the Kolmogorov "K41" scaling to the anomalous (multi-) scaling - 
singular dependence of the Galilean invariant quantities (like e.g. equal-time structure 
functions) on L, characterized by an infinite set of independent anomalous exponents. 
These effects lead to infrared (IR) divergences in the diagrams of perturbation theory 
which, formally speaking, manifest themselves as poles at some finite values of e. Such 
poles and their physical interpretation were discussed in a number of papers; see [7] [TO]. 

The aforementioned crossovers have no analog in the models of critical behaviour, 
but such IR poles are present in their diagrams, too. Moreover, they approach closer and 
closer to the origin e = when the order of the perturbation theory (complexity of the 
corresponding diagrams) increases. The correct resummation of these IR singularities 
is accomplished by the short- distance operator-product expansion (SDE); it shows that 
the finite limit L ^ oo (where L is some IR scale) in the correlation functions exists 
and the singular-in-L terms give only subleading corrections to the scaling behaviour 
[2], [3] . Thus the existence of IR singularities at finite values of e does not hinder the use 
of the RG technique and the e expansion for the description of critical behaviour. 

The SDE technique is equally applied to the stochastic NS problem. The 
distinguishing feature, specific to models of turbulence, is the existence in the 
corresponding SDE of composite fields ( "operators" ) with negative scaling dimensions. 
Such operators were named "dangerous" in [11] because their contributions to the 
OPE diverge for L oo. The summation of the most singular contributions coming 
from the operators f", powers of the velocity field, was performed in Ref. [H], the 
generalization to the case of a time-dependent large-scale field is given in [12]. This 
gives the adequate description of the sweeping effects within the RG formalism; see also 
the general discussion in Refs. [5l [6]. 

According to the SDE scenario, anomalous multiscaling in the structure functions 
can be related to the existence of Galilean invariant dangerous operators. This idea 
was successfully realized for the model of a scalar impurity field passively advected by 
the Gaussian velocity field with given correlation function oc 6{t — t')/k'^~^'^, known as 
the Obukhov-Kraichnan rapid-change model. In the zero-mode approach, which can 
be viewed as a variant of the self-consistency equations, the anomalous exponents were 
derived analytically to order 0(e) [13] and 0(1/ d) [13]; see also [15] for a detailed 
discussion. In the RG and SDE approach, the exponents are identified with the scaling 
dimensions of composite operators built of the scalar gradients (namely, powers of the 
dissipation rate of the scalar field fluctuations, and their tensor analogs) and calculated 
within the expansion in the exponent e (also denoted by ^) to order (three-loop 
approximation of the RG); see the original works [16], [17] and the review paper [18] for 
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generalizations and more references. 

In the stochastic NS problem, however, all the critical dimensions are strictly 
positive for small e and can become negative only for some finite values of e. Therefore, 
dangerous invariant operators cannot be identified within the e expansions and it is 
desirable to construct an alternative perturbation scheme valid for finite e. Attempts 
were made to modify the model by introducing replicas of the velocity field and to 
construct an expansion in [19], but such modifications were inconsistent with the 
Galilean symmetry [20] ; see also discussion [21], [22] for the case of discrete shell models. 

Thus we naturally return to an old idea of the expansion in 1/d, which has 
repeatedly been introduced in various contexts in turbulence [23l [Ml [25l [26l [271 EH] 
and looks very attractive for a few reasons. 

One can hope that in the limit d ^ oo intermittency and anomalous scaling 
disappear or acquire a simple "calculable" form and the finite-dimensional turbulence 
can be studied within the expansion around this "solvable" limit [23]. However, 
in contrast to the 0(iV)-symmetric model of the critical behaviour, where the limit 
^ oo is described in a closed form by the exactly solvable spherical model [29], 
no drastic simplifications were found in [21] for the infinite-dimensional turbulence: all 
classes of diagrams for the velocity correlation function (for nonstationary perturbation 
theory with the Taylor expansion in time and for renormalized perturbation theory 
with skeleton diagrams with dressed lines) survive in the large-c? limit, and the 
incompressibility condition (and thus the nonlocal pressure effects) remain important. 
Much later it was argued (on the basis of a certain SDE-motivated ansatz for dissipative 
terms) that the K41 theory becomes exact and the multiscaling indeed disappears for 
d = oo [25], as it also happens for the Obukhov-Kraichnan model [23]. What is more, 
for the latter it was possible to find the 0{l/d) contribution to the anomalous exponents 
[H]. Numerical simulation of the passive scalar advection in the Obukhov-Kraichnan 
model for large d (up to d = 30) was performed in [30]. However, the systematic 
expansion in l/d has not been yet constructed for that model, let alone the stochastic 
NS equation. 

The key idea of the present paper is to combine the large-d limit with the RG 
approach and the expansion in e. It was noticed earlier in [27l |2H] that taking the limit 
(i — > oo leads to serious simplifications in the RG calculations. In a very important 
paper [27] , scaling dimensions of all the powers of the local energy dissipation rate were 
calculated for d = oo to first order in e. For a finite d, such calculation becomes an 
extremely daunting task due to the mixing of operators in the renormalization, and it 
was not completely performed even for the second power of the dissipation rate; see 
[3T] . The scaling dimensions calculated in [27] appear positive for £ < 2 and vanish at 
the physical value e = 2, in agreement with the arguments of Refs. [211 [25] that the 
K41 scaling exactly holds for ci = oo. It is then possible that the 0{e/d) correction to 
the result of [27] will make the dimensions negative and the anomalous scaling will be 
established as a fact reliable within the double expansion in e and 1/d (we recall that the 
powers of the dissipation rate are the operators "responsible" for the anomalous scaling 
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in the Obukhov-Kraichnan model (THl [13 [IB])- For the vector analog of the rapid- 
change model, where the calculation of the anomalous exponents also becomes rather 
involved already in the 0{e) order due to the mixing of composite fields, additional 
expanding in 1/d also leads to serious simplification ^28j . 

In this paper we systematically investigate the infinite-ci limit within the RG 
framework. Although we were not able to exactly solve the problem for d = oo nor 
to construct an analog of the "spherical model," we obtained a number of new results 
which we consider interesting and which, as we believe, can give some hints for the 
further study of the large-d behaviour. 

In section [2] we briefly recall the field theoretic formulation, renormalization and 
the main findings of the RG approach for the stirred NS equation. A more detailed 
specification of the limit (i — cxo is given. In section [3] we consider the linear response 
(Green) function and show that, for d oo, lots of diagrams vanish and the remaining 
ones simplify drastically. This allows us to analytically calculate in the three-loop 
approximation the key ingredients of the RG approach: the renormalization constant, 
f5 function, coordinate of the fixed point and the ultraviolet (UV) correction exponent 
uj (for finite d, only two-loop results for these quantities are known, and they involve 
numerical calculations of the corresponding diagrams ^Q\)- Some speculations are also 
made about the hypothetical expressions for these quantities beyond the scope of the e 
expansion. 

The two-point velocity-velocity correlation function is analyzed in section |H Here, 
no remarkable simplifications in the diagrams occur in the limit d — >^ oo, in agreement 
with the observations made earlier in Ref. [21]. However, we propose some trick which 
allows us to perform the analytical calculation of the pair correlator in the two-loop 
approximation (which, in orders of the coupling constant, corresponds to the three- 
loop approximation for the Green function). Although the results can be interpreted a 
posteriori in terms of the steepest-descent method, the practical use of the latter beyond 
the one-loop approximation is hardly possible due to complexity of the corresponding 
integrals. 

These results are exploited in section [51 There, we calculate to the third order in e 
the Kolmogorov constant Ck in the spectrum of turbulent energy and the inertial-range 
skewness factor S. An important point here is not only the inclusion of the third-order 
correction, but also the derivation of Ck through a universal (in the sense of the theory 
of critical behaviour) quantity. This approach was proposed earlier in [10|, where Ck 
was calculated (for d = 3) to the second order in e. It allows one to obviate the main 
shortcoming of earlier calculations of Ck'- the intrinsic ambiguities in the corresponding 
e expansions (see e.g. [32] and the discussion in [6]). The results are briefly summarized 
in the Conclusion. 
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2. Stochastic Navier— Stokes equation, choice of the random force and the 
field theoretic formulation 

As the microscopic dynamical model of the fully developed, homogeneous, isotropic 
turbulence of an incompressible viscous fluid one usually takes the stochastic NS 
equation with a random stirring force 

VtVi = uod\~diV + fi, Vt = dt + vA, (2.1) 

where Vi is the transverse (divergence-free, due to the incompressibility condition 
diVi = 0) velocity fleld, V and fi are the pressure and the transverse random force 
per unit mass (all these quantities depend on x = {t, x}), uq is the kinematic viscosity 
coefficient, d"^ is the Laplacian and Vt is the Lagrangian derivative. The problem (12.11) 
is studied on the entire t axis and is augmented by the retardation condition and the 
condition that Vi vanishes for t — oo. 

We assume for / a Gaussian distribution with zero mean and correlation function 

{fiix)f,ix')) = J dkP,,{k)df{k) exp[ik(x-x')], (2.2) 

where Pij(k) = 6ij — kikj/k^ is the transverse projector, df{k) is some function of the 
wave number A; = |k| and model parameters, and d is the dimension of the x space. The 
time decorrelation of the random force guarantees Galilean invariance of the stochastic 
problem ^J^, flO) . 

Let us specify the form of the function df{k) in the correlator fl2.2l) used in the 
RG theory of turbulence. Physically, the random force models the injection of energy 
into the system due to interaction with the large-scale motions. Idealized injection by 
infinitely large vortices corresponds to 

df{k) = 2(271^8 6(k)/{d-l) , (2.3) 

where S is the average power of the injection (equal to the average dissipation rate) and 
the amplitude factor comes from the exact relation 

^ = 7^ / d^dfik). (2.4) 



2(27r)"' 

On the other hand, for the use of the standard RG technique it is important that the 
function df{k) have a power-law behaviour at large k. This condition is satisfied if df{k) 
is chosen in the form 

df{k) = Dok^-'^-^', (2.5) 

where Z^o > is the amplitude factor and £ > is the exponent with the physical value 
e = 2. This can be explained with the aid of the well-known power-law representation 
of the (i- dimensional S function 

6{k) = hm ^ I rfx(Ax)^-^ exp[i(kx)] = 

= S^'k-'^hm[{4~2e){k/Ar-'^], (2.6) 
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with some UV momentum scale A. Here and below we denote 

S, = 27r'^/Vr(d/2), Sd = Sd/{27rf, (2.7) 

where Sd is the surface area of the unit sphere in c?-dimensional space and r(. . .) is the 
Euler Gamma function. For e ^ 2, the function (12. 5p turns to the ideal injection (12.31) 
if the amplitude Do is related to S as 

_^ 4(2_- e)A^'-^ ^ ^^^^^ 2. (2.8) 

As already mentioned, in the RG approach to the problem (12.11) . (12.21) . (12.51) the 
exponent e plays the part analogous to that played by 4 — d in Wilson's theory of 
critical phenomena, while d remains a free parameter. A more realistic model is 

dj{k) = Do k'^-'^-^' h{m/k), h{0) = 1, (2.9) 

where m = 1/L is the reciprocal of the integral turbulence scale L and h{m/k) is 
some well-behaved function that provides the IR regularization. Its specific form is 
unessential; we shall always use the sharp cutoff 

h{m/k) = e{k-m), (2.10) 

which is the most convenient choice from the calculational viewpoints; 6(. . .) being the 
Heaviside step function. 

According to the general theorem [33], stochastic problem (12.11) . (12.21) is equivalent 
to the field theoretic model of the doubled set of fields $ = {v', v} with action functional 

S{<l>) = v'D,v'/2 + v' [-Vt + i^od^] V, (2.11) 

where D^ is the correlation function (12. 2p of the random force fi with the function df{k) 
from (12. 9p and all the required integrations over x = {t, x} and summations over the 
vector indices are understood, for example. 



v'{vd)v = J dt J dx vl{vjdj)vi. 

The auxiliary vector field is also transverse, 9jf ■ = 0, which allows to omit the pressure 
term on the right-hand side of Eq. (12.111) . as becomes evident after the integration by 
parts: 

Of course, this does not mean that the pressure contribution can simply be neglected: the 
field v' acts as the transverse projector and selects the transverse part of the expressions 
to which it is contracted in (12. lip . 

Formulation ( 12.111) means that statistical averages of random quantities in the 
original stochastic problem (12.10 . (12.20 can be represented as functional averages with the 
weight exp5'($), and the generating functionals of total [G'(yl)] and connected [Vr(A)] 
correlation functions of the problem are represented by the functional integral 

G{A) = exp W{A) = J'^'^ exp[5($) -|- A$] (2.12) 
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with arbitrary sources A = {A^' , A"} in the hnear form A$ = J dx A^{x)^{x). 

The model (12.111) corresponds to a standard Feynman diagrammatic technique; the 
bare propagators (hues in the diagrams) in the time-momentum (t-k) representation 
have the forms 

{v,{tmt'))^ = 0(t - t') exp {-lyo^it - t')} P,,(k) = h, 



(^.(^)^.(^'))o = ^ exp {-uok'\t - t'\} P,,(k) = , 

(^:W(0)o = o/ (2.13) 

with df{k) from l\2.9h and the Heaviside step function 6(. . .). The interaction in (I2.1ip 
corresponds to the triple vertex —v'{vd)v = with vertex factor 

Vl = iik,6,s + kAj) = (2.14) 

where k is the momentum argument of the field v'. 

It is convenient to introduce the new parameter ("coupling constant") by the 
relation 

90 = Do/ul (2.15) 

so that go oc A^^ with the UV momentum scale A coming from (12.61) . Thus the model 
(12. lip becomes logarithmic (the coupling constant becomes dimensionless) at e = 0, and 
the UV divergences manifest themselves as the poles in e in the correlation functions 
of the fields $ = {v,v'}. Dimensional analysis (power counting) shows that superficial 
UV divergences, whose removal requires counterterms, can be present only in the 1- 
irreducible correlation functions {v'v)i-i^ and the corresponding counterterms reduce to 
the forms v'd'^v, v'dfV and v'{vd)v. In fact, in the model (12.1 ip there are fewer possible 
counterterms than allowed by the naive power counting. The symbol d at the vertex 
in (12. lip can be moved onto the field v' using the integration by parts, which means 
that the counterterms to the 1-irreducible functions must contain at least one spatial 
derivative per each field v'. This excludes the structure v'dfV. The Galilean symmetry 
requires that the counterterms v'dtv and v'{vd)v must form the invariant combination 
v'VtV, which excludes the structure v'{vd)v. In the special case ci = 2 a new UV 
divergence appears in the 1-irreducible function {v'v')i-ij.. Since we are interested in the 
large-d limit here, we can simply ignore that divergence. 

Then the inclusion of the only remaining counterterm v'd'^v in the action functional 
(12. lip is reproduced by the multiplicative renormalization of the parameters uq and go 
with the only independent renormalization constant Z,^: 

1^0 = uZ,, go = gfi''Zg, Zg = Z;' {Do = gou^^ = af^'^'^')- (2.16) 

Here /i is the reference mass in the minimal subtraction (MS) scheme, which we always 
use in what follows, g and u are renormalized analogs of the bare parameters go and uo, 
and Z = Z{g,6,d) are the renormalization constants. No renormalization of the fields 
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and the IR scale niQ = m is needed, i.e., Zq, = 1 for all $ and = 1. The renormalized 
action has the form 

5(<l>) = v'D,v'/2 + v' [-Vt + i^Z^d^] V, (2.17) 

where the amplitude Dq in Dy is expressed in renormalized parameters using the last 
relation from fl2.16p . 

In the MS scheme the renormalization constants have the form "1 + only poles in 
e" in particular, 

oo oo n 

Z, = l + Y, ak{g)e-' = 1 + 5^ ^7" 5Z «nfc£"', (2.18) 

k=l 71=1 k=l 

where the coefficients a„fc depend only on d. 

Since the fields are not renormalized, their renormalized correlation functions 
coincide with their unrenormalized analogs G = ($...$); the only difference is in 
the choice of variables and in the form of perturbation theory (in g instead of ^fo)- 
G^((7, 1/, /i, m, . . .) = G{go,h'o,mo, . . .). Here the dots stand for other arguments like 
coordinates, times, momenta and so on. We use T>^ to denote the differential operator 
fid^ for fixed bare parameters go,i'o, niQ and operate on both sides of that relation with 
it. This gives the basic differential RG equation: 

'DRcG^'ig, u, /i, m, . . .) = 0, Vrg = 'D^ + (3{g)dg - i,{g)V,, (2.19) 

where 'Drg is the operation T)^ expressed in renormalized variables, Vx = xdx for any 
variable x, and the RG functions (the anomalous dimension and the (3 function) are 
defined as 

7^(5f) = In Zy = P{g, e)dg In Z^, 

Pig, e) =V^g = g [-2e + 87.(^7)] ; (2.20) 

the relation between (3 and 71, results from the definitions and the last relation in fl2.16p . 
Combining the two relations in fl2.20p and substituting (12.181) gives 

lvi.9) = ~'^1^gO'i{g) + the terms containing only poles in e. (2-21) 

with ai{g) from (12.181) . From the UV finiteness of the renormalized functions G^ it 
follows that the pole terms in (I2.2ip cancel each other, which eventually gives 

00 

7,(^7) = -2Vga,{g) = -2 na^if?" (2.22) 

n=l 

with ai{g) and a„i defined in (12.180 . The cancellation occurs due to the fact that the 
coefficients ank with k > 1 are algebraically related to a„i, the coefficients in front of 
the first-order poles, for example, 022 = a^i- Such relations can be used to check the 
results of the practical calculations. 
The one-loop result 

an = -(rf- l)Sd/8(rf + 2) (2.23) 
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with Sd from (12.71) is well known, while the two-loop coefficient 021 was calculated (for 
a few values of d, including d = 3 and in the limits d ^ 2 and d 00) in [10]. In the 
latter limit, which we are interested in here, one obtains: 

a2i = a?i |-^ + 0(l/rf)|. (2.24) 

From the ffist order result (I2.23P it follows that the (3 function has a nontrivial fixed 
point [(3{g*) = 0] in the physical region g > with the coordinate 

g, = 8{d + 2)e / 3{d - 1)3^ + 0{e^). (2.25) 

The correction exponent u = [3\g^) = 2e + Oi^e^) > at this point is positive, so that 
it is IR attractive and governs the IR behaviour of the correlation functions. The value 
of 7iy(5') at the fixed point is found exactly using the relations (12.201) : 

ll = iM = 26/3. (2.26) 

The relations (12.231) and (12.241) illustrate the fact that the coefficients that have a 
well-defined finite limit for d ^ 00 are bni = ctni "S*^" rather than a^i themselves (this 
is obvious from the practical calculation, see the next section [3]). It is convenient to 
introduce the new coupling constant by the relation 

u = gSd. (2.27) 

Then the correlation functions G and G^, renormalization constants (I2.18P and the RG 
functions 7,^ and (3 = T>^u will have well-defined finite limits d ^ 00 when expressed as 
perturbation series in the parameter u; the latter is kept finite in that limit. 

In the following, the large-c? limit will always be understood in this sense. From the 
physics viewpoints it corresponds to the natural requirement that the quantity which 
is kept finite for — cxd is the energy input per one component of the velocity field 
(that is, per one spatial dimension), which is clear from the relations (12. 7p . (12. 8p and 
( 12.15p . Such a choice is also consistent with further RG analysis which shows that the 
fixed-point value u^, of the parameter u (and not that of g) is finite for d — 00, as 
illustrated by equation (12.250 . 



3. Calculation of the renormalization constant and RG functions to the 
third order 

Let us turn to the calculation of the renormalization constant Z,^ in (12.180 with the 
accuracy 0{g^) (three-loop approximation) in the limit d — 00. In general, the 
renormalization constant Zr corresponding to a certain 1-irreducible Green function 
can be found from the relation [3l] (see also the monograph [3]) 

Z^=l-/C7^'f, (3.1) 

where F is the function F normalized with respect to unity in the zeroth order of the 
perturbation theory, TZ' is the incomplete 7^-operation which involves all the subtractions 
for the divergent subgraphs, without the last subtraction for the diagram as a whole. 
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and /C is the subtraction operation for the given renormahzation scheme. In the MS 
scheme, which we use in our calculation, /C subtracts only poles in e: 

oo —1 
n=—oo n=—oo 

for any Laurent series. The function F in (13. ip should be calculated in terms 
of the renormalized parameters from the "basic" action functional [SI [35] which is 
obtained from the renormalized action (12.171) by the replacements Zi —* 1 for all the 
renormalization constants. 

The only independent renormalization constant Zj, in our model (12.171) is 
determined by the 1-irreducible function (in the frequency-momentum representation) 

Tij{uj, p) = {v[vj)i^i, = T{uj, p)Pij(p), (3.2) 

where the projector -Pij(p) = Sij — PiPj/p"^ arises due to the transversahty of the fields 
f Vj. Thus the scalar coefficient r(ci;, p) in (13.21) is given by the relation 

T{uj,p) = r,,(cu,p)P,,(p)/(rf- 1) = Tu{uj,p)/{d-l), (3.3) 

where {d—1) = -Pm(p) comes from the trace of the projector. In the perturbation theory, 

T{uj,p) = iuj - up^ + 0{g). 

Since the counterterm to the function (13.21) is proportional to vp^ , one can put = 
in r and neglect the higher-order terms in p^. Therefore the normalized scalar function 
r in ( 13. ip can be taken as 

r=lim ^ = 1 + 0^; 3.4 

it depends only on the completely dimensionless variables g and m//i. 

The elementary integrations over times (or equivalently frequencies) in the diagrams 
are always easily performed; the serious problem is only the integration over the 
momenta. (The factors u arising from these integrations and from the propagator {vv)o 
group together to cancel u in the denominator of ( 13. 4p .) 

The first important observation is that, for any diagram of the function (13. 2p . the 
number of loops (that is, the number of independent integration momenta) is equal to 
the number of the lines (propagators) of the type {vv)q (this is easy to understand from 
a few examples, given in this section, so that we shall not give the formal proof). It 
is therefore possible to assign the "pure" independent momenta (which we shall denote 
by k, q and 1 for the three-loop diagrams) to these propagators. Then the momenta 
associated with the remaining lines of the type {vv')q will be certain definite linear 
combinations of the external momentum argument p and the integration momenta k, 
q etc. 

Each line {vv)o brings about one integration over an independent pure momentum: 
;{vv)o{k) ... = / dke-''-^' .... (3.5) 



{27tY'"'' (27r) 
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Here we used the last expression from fl2.13l) and omit the dimensionless time-dependent 
factor and the projector, denoted by the elhpsis. Now we can separate the integration 
over k into the integrations over the modulus k = |k| and the direction n = k/fc 
("angles"), which gives dk = k'^~^dkdn, and replace the integration over the angles by 
the averaging {■ ■ ■)s over the unit (i-dimensional sphere, normalized such that (l)^ = 1. 
This gives 

1^ dk k'~'^ j dn... = g^i'^u'S, P dk k'-'^ {. . .)s. (3.6) 

The factor Sd from (12.71) groups together with the coupling constant g; this quantity is 
kept fixed in the limit d oo. The lower limit m in the integration over k comes from 
(I2.10p and provides the IR regularization. 

In (13. 6p the powers A;"'^ from the correlator (12.91) and k"^ from the Jacobian cancel 
each other, so that the dependence on d in the exponent disappears. The remaining 
dependence on d comes from the contractions of the projectors (12.131) and vertices (I2.14p . 
with further integrations over the angles. 

The vertex V^^^ in (I2.14p is explicitly symmetric with respect to the indices js, that 
is, symmetric with respect to the two attached fields vj and Vg- Let us split it into the 
two asymmetric parts 

V-^ = -I- i'\ 1\ 

ijs ijs ' * isji ' J 

where the new asymmetric vertex V^^^ = ikj6is in the diagrams will be denoted as 

i IT 

V^^ = ikj6,s = -ipjSis = s^^„^t|;(^*p), (3.8) 
so that the splitting (13. 7p is diagrammatically represented as 

Here, the end marked with the slash corresponds to the field v', the end with no 
marks corresponds to the field v which stands under the derivative in the expression 
—v'{vd)v = v[Vj^gVjVs] the vector indices of these two fields are contracted to each 
other. The end marked with the thick dot corresponds to the second field v (without 
a derivative); its vector index is contracted with the index of the momentum. The 
momentum k is the argument of the field v' . In the asymmetric vertex, it can be 
equally replaced with — p, the argument (up to the minus sign) of the unmarked field 
f , due to the transversality of the marked field v. 

The use of the asymmetric vertex leads, of course, to increase in the number of 
diagrams: each diagram with symmetric vertices (I2.14p gives rise to the sum of 2^ 
contributions, represented by diagrams with vertices (13. 8p . as illustrated in figure 1 
for the one-loop case. This inconvenience is compensated by the fact that the large-o? 
behaviour of such individual contributions is easier to analyze. We shall see below that 
most of them vanish in the limit (i — oo, and it is possible not to draw the corresponding 
diagrams from the very beginning. 
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Figure 1. Passage to the asymmetric vertex in the one-loop approximation. 



From (13.31) it follows that a diagram can survive in the limit d ^ oo only if its 
contribution to the quantity Tii{u!,p) in the numerator behaves as 0{d) for large d. 
The only possible source of such a contribution is a closed (self-contracted) chain of 
intermittent (5-symbols, 



coming from the vertices (13.81) and from the projectors in propagators (I2.13p . with the 
final contraction with 6ji from the projector in (13. 3p . It is easily seen that, in terms of 
the asymmetric vertex, each diagram can involve no more than only one chain of the 
type (13.10p . Assume for definiteness that the factor 6ii-^ in ( 13.10p comes from the leftmost 
vertex of the diagram, then 6i^j comes from the rightmost vertex. Then from the form 
of the vertex (13. 8p it follows that the vector index i corresponds to the field v with a 
derivative, and ii corresponds to the field v'. The latter is necessarily contracted with 
the field v in the neighboring vertex, because the propagator {v'v')q in (12.131) vanishes 
identically. Thus the next factor 6ij^i^ comes from the propagator {v'v)o. It is contracted 
with the next symbol 6i^i^ from the vertex, in which ^2 must correspond to the field v 
with a derivative (any contraction with the momentum in (13.80 would break the chain 
of 5 symbols); then i^ necessarily corresponds to the field v', and so one. 

We may conclude that the product of the type (13.101) is necessarily associated with 
the chain of propagators {v'v)o which begins in the leftmost vertex of the diagram and 
ends up in the rightmost vertex; we shall call it "the spine" in the following||| It can 
be formally proved, but is easily seen from the examples, that each diagram of the 1- 
irreducible function (13.21) involves one and only one spine. However, the spine gives rise 
to the factor of the type (13.101) only if all the fields v enter it with derivatives, so that 
their vector indices are contracted with the indices of the fields v'; this means that all 
the dots in the vertices (13.80 must be placed out of the spine. This fixes the placements 
of the dots in the vertices that belong to the spine in the unique way and drastically 
reduces the number of relevant diagrams. In figure 1, only the first diagram contains 
the factor (13.101) and gives a nonvanishing contribution to the quantity (13. 3p : the other 
diagrams vanish for d ^ oo and can be dropped from the very beginning. Figure 2 gives 

X The factor of the form (|3.10p can also be associated with a closed circuit of the propagators {v'v)o, 
not connected to the external vertices of the diagram. But such diagrams vanish identically owing to 
the fact the propagator {v'v)o is retarded: it contains the step function in the time arguments; see 




(3.10) 



(I2l3ll . 
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Figure 2. Examples of three-loop diagrams which do not vanish for d ^ oo. 



examples of nonvanishing three-loop diagrams with the proper placement of the dots on 
the spine (here and in the figures below, the spine is always shown as a horizontal chain 
of propagators). 

The next important observation is that any nontrivial angular integration brings 
about smallness in 1/d; the corresponding contributions in the diagrams can be dropped 
for d ^ oo. More precisely, the average over the unit sphere of the scalar product of two 
independent momenta k and q with the angle i} between them is given by the relations 



and ((kq)^"^^)^ = 0; here n = 1,2 and so on. In the three-loop diagrams, averages of 
the form 



arise in the integrands. The corresponding expressions, analogous to (13. lip , are rather 
complicated; they can be found e.g. in sec. VIC of Ref. [IT]. It is only important for us 
here, that they always vanish for c? — oo (with the obvious trivial exception = for 
all i.) 

In all the relevant diagrams, the momentum argument of the two external vertices 
can be chosen equal to p, the external momentum argument of the function (13.21) . [This 
would be impossible for the diagrams in which the dot in the rightmost external vertex 
is placed on the external line, like in the last two diagrams in figure 1. The momentum 
arguments for such vertices necessarily involve internal integration momenta k, q and 
so one. But we have already established that for d ^ oo their contributions vanish.] 
Thus the overall factor of the form piPj is isolated from the quantity (13.31) . and, since 
we are interested in the limit (13. 2p . we can set p = in the rest of the integrand. 

If one of the external momenta, say Pi, is contracted with some integration 
momentum in the integrand, say ki or qi + ki, a construction with scalar products 



((kq)2")5 = (fcg)2"(cos2"^9)s = 




0(rf-") 



(3.11) 



{{kqrikiniqDs 



(3.12) 
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Figure 3. Examples of three-loop diagrams which have no "second spine" and vanish 
for d —^ oo. 



of the form fl3.1ip or fl3.12p arises, and in the hmit d — > oo the diagram vanishes after 
the angular integrations. Nonvanishing contribution arises only if the external momenta 
are contracted with each other through a chain of S symbols: 



In the diagram, it corresponds to the chain of propagators that connects the two 
external vertices, with only two dots placed on the chain at the external vertices; they 
correspond to the external momentum. Each dot placed on the chain at an internal 
vertex corresponds to contraction with some integration momentum k, q etc, and 
therefore breaks the chain fl3.13p . All the diagrams shown in figure 2 have such "second 
spine," having the form of a big arc connecting the external vertices and surrounding 
the diagram from above; the proper placement of the dots is also shown. From these 
examples it is clear that the second spine always has the form 



that is, it consists of two chains of mixed propagators (ff')o and {v'v)o (not necessarily 
of equal "lengths"), proceeding from the opposite directions and connecting with the 
only propagator {vv)q. In the simplest case of the only nonvanishing one-loop diagram 
in figure 1, the "second spine" exists and consists of the single propagator {vv)q. 

If such "second spine" exists, it is unique because the chain of 5 symbols in (13. 130 
"cannot branch." But the existence is not guaranteed (in contrast to the "first spine" 
discussed above), as illustrated by the diagrams in figure 3. The first spines, shown as 
horizontal lines, consist of 5 and 3 propagators, respectively (the proper placement of 
the dots, which gives rise to the chains of the form fl3.10p . is always unique and is not 
shown). But the second spines do not exist: for the first diagram, there is no chain 
of propagators, connecting the external vertices. For the second diagram, such a chain 
exists (the big arc) but it involves two propagators {vv)o and thus it does not have the 
form f l3.14p : any placement of the dot in the right upper vertex breaks the chain of 5 
symbols in fl3.13p . 

Such diagrams give no contribution for — * oo, so they can be dropped from 
the very beginning. In the nonvanishing diagrams, the factor coming from the chain 
fl3.13p . explicitly cancels the analogous factor in the denominator of (13.40 : the remaining 
quantity depends only on the dimensionless parameters g, e and ii/m (the latter will 
also disappear soon). 




(3.13) 




(3.14) 
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Discarding the diagrams which vanish for d oo (those with "wrong" placement 
of the dots on the spine or those having no second spine) leads to impressive decrease in 
the number of the surviving diagrams. In the one-loop approximation only one of four 
diagrams survives (as is clear from figure 1), in the two-loop approximation 6 diagrams 
of 120 survive, and in the three-loop approximation only 80 diagrams survive of 8160. 

The final simplification in the remaining nonvanishing diagrams is the possibility 
to drop all the scalar products of the integration momenta k, q and so on, due to the 
smallness in the averages in expressions (13.111) and (13.121) . If a scalar product appears 
in the denominator of the integrand, say k"^ + q"^ + (kq), these expressions are not 
applied directly, but one can expand all the denominators in all the scalar products 
and then discard all the terms except for the first one (with no scalar products) due 
to relation ( 13. lip and its analog for more complicated structures like ( 13.12^ . This is 
the same as to drop all the scalar products in the integrals from the very beginning. 
The diagrams in which the integrands are proportional to some scalar products vanish 
after this procedure, but we have no general rule how to identify them from the very 
beginning. The averaging over the angles now becomes trivial: (1)5 = 1. 

Thus the expressions for the remaining nonvanishing diagrams in the n-loop 

approximation are represented by n-fold integrals over the n scalar variables, the moduli 

of the independent integration momenta k, q and so one. It is convenient to pass to new 

integration variables fc^, and so one; then will appear as the new IR cutoff in the 

integrals. We shall retain the same notation k, q, m for these new parameters. Then, 

for n = 3, such integrals take on the forms 

dk dq dl 

J{m)= -j^ f^k,qj), (3.15) 

J m Jm H J m 

where / is some dimensionless function. From dimensionality considerations it is clear 
that the integral (13.151) is the definite power of the IR scale, J{m) = m~^^J(l). Then 
applying the differential operation = md/dm to the integral (I3.15P multiplies it by 
the factor {—3e), and the following identity holds: 

J(l) = -^P^J(m)U=i. (3.16) 

The parameter m enters the expression ( 13.15^ only in the lower integration limits, so that 
explicit differentiation in the right-hand side of (13.160 reduces the number of integrations 
by one. Since one pole in e in (I3.16P is already isolated explicitly, in the resulting double 
integrals one should retain the terms of order e~^, and e". 

The complete three-loop calculation involves only six basic integrals Jj of the form 
()3.15p with the following functions ff. 

n _ kql _ kq kql 

h — 71— , — ri7\-> h — 71—, — riv?^ 



(k + q + l)^' {k + q + iy {k + q){k + l){q + iy 

^ kq ^ kq ^ kq 

{k + q){k + q + iy {k + q){k + iy {k + qf' ^ ' ' 



The practical calculation gives 



-1 -1 ~ _ 7r2 -1 

— — , Jo — , 0% — , Ja — ~l~ 1, 

2' 2e 4 ' 2£ ' 
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J5 = ^ + y, J6 = |-31n2, (3.18) 

where we denoted Jj(l) = (l/Se) Jj. 

The 7^-operation in fl3.1l) imphes subtraction of the UV divergences related to the 
subgraphs; this is a standard routine, see e.g. [Sj [35]. It requires calculation of the 
one-loop and two-loop analogs of the integrals like (13.151) with the proper accuracy; this 
is of course an easier task than the calculation of the three-loop integrals. Combining 
all the contributions finally gives the following third-order result for the renormalization 
constant (l238l) : 

Z -1 ^ ^'(^ + ^) n3(10 + 9. + 7.^) 

with the coupling constant u introduced in (I2.27p . It is worth noting that the 
transcendental numbers (In 2 and vr^) entering the expressions (13.181) have disappeared 
in (I3.19P as a result of the subtraction of sub divergences. 

Then for the anomalous dimension (12.201) using the relation (I2.2ip one finds: 

Mu) = ^u + ^u' + ^u' + Oiu') (3.20) 

and the three-loop (3 function is found by the second relation in (I2.20p . Solving the 
equation /3(m*) = perturbatively in e, for the coordinate of the fixed point one obtains: 

u. = le-le'-P' + 0{e'). (3.21) 
The UV correction exponent 

9 1 

u; = p'{u,) =2e + -e' + -^e' + 0{e^) (3.22) 

completes the hst of the third-order results. 

The expressions obtained appear rather simple: in particular, all the coefficients 
in the series in or e for Z^, (3, and uo are rational numbers, although for 
dynamical models such coefficients usually involve transcendental quantities (logarithms, 
hypergeometric functions etc) already in the two-loop approximations. One can hope 
that this fact remains valid for all the higher-order terms. It is then tempting to try to 
guess their general form and, therefore, to construct some hypothetical exact (that is, 
beyond the expansion in g or e) expressions for those quantities. Such exact answers 
for all the quantities (I3.19p - (l3.22p exist |36] for the well-known Heisenberg model of the 
developed turbulence |37|, |38j. In particular, for uj and (in the MS scheme) one has 
BSl: 



cu = £(l-£/3)/(l- 5^/12), -u, = (e/3)(l -£/3)^/^ (3.23) 

where all the parameters have the same meaning as in our model (it is implied, of course, 
that the usual Heisenberg model is generalized to the energy pumping of the form (12. 5p ). 
Let us assume that in our model the fixed-point coordinate has a similar form: 

= {Cie){l +£02)'"'. 
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Then comparison with the cut (13.211) of the e expansion for u^, allows one to find all 
the coefficients q, which appear very simple and, surprisingly enough, the exponent 
C3 = 1/4 is exactly the same as in the Heisenberg model: 

= (8e/3)(l - Ae/3y/\ (3.24) 

Now the exact expression for u can be find from (I3.24p . Differentiating the equation 
(I2.20p for the function /3 = "D^u with respect to u, setting u = u^, and taking into 
account 7^ gives: 

u = (3'{u^) = 37^(u*), 

and differentiating the equation (12.261) with respect to e gives: 

llM<{e) = 2/3 

(the first differentiation with respect to u and the second - with respect to e). Combining 
these equalities gives the exact relation 

uj{e) = 2uju'^. 

Substituting the explicit expression (I3.24p for the function u^: = «*(£) gives the desired 
hypothetical exact answer for the UV correction exponent: 

= 2^(1 -4e/3)/(l -5£/3), (3.25) 

which looks very similar to the corresponding expression (I3.23P for Heisenberg's model 
and coincides (as is easily checked) with the expansion (I3.22p . From (I3.24p and (I3.25P it 
is then possible to derive the exact expressions for the (3 function, anomalous dimension 
and the renormalization constant Z^; cf. [36] for the Heisenberg model and [39] 
for the Kraichnan model. These expressions are rather cumbersome and will not be 
reported here. 

Expressions (13.240 . (13.250 . provided they are correct, show that the IR fixed point 
> disappears at e = 3/4, that is, before the physical value e = 2 is achieved. This 
must lead to some qualitative change in the IR behaviour of the correlation functions 
which cannot be properly described within any finite-order approximation of the e 
expansion. Of course, one should not take too seriously these hypothetical expressions, 
which, of course, are not the only possible ones. 

4. Pair correlation function: Third-order approximation 

In this section we shall discuss the equal-time pair correlation function of the velocity 
field in the momentum representation, 

Dij{p) = P,,{p)D{p). (4.1) 

We are ultimately interested in its inertial-range behaviour (m ^ p), so in the following 
we set m = 0; the IR regularization in the diagrams is provided by the momentum p. 
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Solving the RG equation (12.191) for the function D{p) gives (see e.g. [21 El E] for the 
detailed derivation) 

D{p) = gu'p-'+'R{s, g) ^ gl'^p-'^'^^ R{1, g,), s = p/fi. (4.2) 

The first equality, along with dimensionality considerations, introduces the function R 
which depends on the two dimensionless variables g and s = p/ fi (dependence on e and 
d is always implied). The second relation holds in the IR asymptotic region s = p/ fi <^ 1 
and involves the coordinate of the IR attractive fixed point g^, and the critical dimension 
of the velocity field, which is known exactly: A^, = 1 — 2e/3. It depends only on the 
amplitude Dq from (12.51) and is independent of the viscosity coefficient uq in agreement 
with the second Kolmogorov hypothesis; see the discussion of this issue in [H [5l El [TT]. 

The function R can be directly calculated in the renormalized perturbation theory 
as a series in g, its coefficients being finite aX e —>■ 0. [One factor g in (14. 2 p is explicitly 
isolated from R such that its expansion in g begins with 0{g^).] Substituting as 
the series in e (with the first-order term given by (I2.25P ) gives the e expansion of the 
amplitude R{l,g^). It was calculated earlier in the one-loop approximation for certain 
values of d; see [10]. Here we shall calculate it in the limit ci — cxd in the two- loop 
approximation (three terms of the e expansion). This accuracy is consistent with the 
three-loop calculation of the renormalization constants and the coordinate g^, ~ 
presented in section El 

The renormalized perturbation theory for the function (14. ip has the form 

with some dimensionless coefficients q = Ci{e,d). The first-order term is the bare 
propagator {vi{t)vj{t'))o from (I2.13p . in which we set t = t' and replaced z/q — > z/ and so 
on. The contribution D^"'^ of a certain ra-loop diagram to the series (14. 3p has the form 

^'"'(P) = (f ) . (4.4) 

the problem is to calculate the dimensionless coefficient A^^^ = A^'^\e,d). After the 
simple integrations over the time variables, the quantity (14. 4p is represented as the 
integral over n independent momenta. 

In contrast to the Green function (13. 2p . now it is impossible, in general, to assign 
these momenta to all the propagators {vv)q. The problem arises if the diagram includes 
a 1-irreducible subgraph of the type (f'f')i_ir, which itself is not a subgraph of another 
1-irreducible subgraph (then it depends as a whole on the external momentum p). Then, 
as illustrated by figure 4, some of the propagators {vv)o necessarily correspond to certain 
linear combinations of the independent momenta and the external momentum p, and 
not to a "pure" independent momentum. Thus the cancellation of the powers like k~'^ 
from the correlator (12. 9p and k"^ from the Jacobian (discussed for the Green function 
below the equation (13. 6p ) does not occur, the dependence on d in the exponent does 
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Figure 4. Examples of the one-loop and two-loop diagrams for the equal-time pair 
correlation function. 

not disappear, and the calculational techniques for c? — oo, developed in the previous 
section, are not directly applicable. 

The general situation is illustrated by specific examples given in figure 4. The 
external momentum is denoted by p and the independent integration momenta are k, 
q and qi,2 (the difference between "k-momenta" and "q-momenta" will become clear a 
bit later). The problem of the non-cancellation does not occur for the diagrams a and 
c: the diagram a has no subgraphs of the type (f'f')i-ir, while the diagram c has one 
such subgraph, but it belongs to another subgraph. Thus it is possible to assign the 
pure integration momenta (denoted by q) or the external momentum p to all the {vv)q 
lines, as shown in the figure. The diagrams b, d and e involve one "harmful" subgraph 
of the form (u'w')i_ir, and one of the {vv)q lines necessarily corresponds to certain linear 
combination of "pure" momenta; in all the examples, it is p — k. These examples also 
illustrate the general fact that any diagram can involve no more than one "harmful" 
{v'v')i_ij. subgraph. What is more, as illustrated by the diagram e, the problem of the 
non-cancellation occurs only within this subgraph; the other subgraphs involve only 
momenta of the type q. 

In principle, the problem can be attacked by the steepest descent method, but in 
practice it is hardly possible to find the stationarity manifold in the integration region. 
In order to circumvent this difficulty, we will apply the following trick. 

Assume that the diagram indeed contains one "harmful" subgraph of the type 
(f'f')i_ii.. Since the p-dependence of the quantity D^'^\p) in (14. 4p is known, we can 
integrate the both sides of that equation over p with an arbitrary "weight function" 
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F{p'^) (the only requirement is that the integral be convergent) and then extract the 
amplitude A^"-^ from the identity 

^ ^ ~ Saj^dpp^-^^(-+')F{p^y ^ ""^ 

In the denominator we have performed the angular integrations, which gives the factor 
Sd from (12 ■71) . The numerator is represented by the integral over {n + 1) momenta: p 
and n independent momenta from the original diagram. As is clear from the examples 
in figure 4, the number of the lines {vv)o in any n-loop diagram is equal to (n + 1). 
Thus it is possible to perform the change of variables in the integral such that the 
pure independent integration momenta will be assigned to all the {vv)o propagators. 
Now the powers of the integration momenta k~'^ from the correlators (12. 9p and k'^ from 
the Jacobians cancel each other, the remaining dependence on d comes only from the 
contractions of projectors and vertices with further angular integrations. Thus we arrive 
at the integrals of the type discussed in the previous section, and they can be calculated 
for (i — > oo in the same manner. In particular, all the scalar products can be dropped, 
and the angular integrals in the numerator of (14. 5 p give only the factor S^~^^; we are left 
with an (n + l)-fold integral over the moduli. 

Let us derive a more explicit expression for the resulting integrals and check their 
independence of the weight function F{p'^). In the new variables the former external 
momentum is certain linear combination of the part of the new independent momenta: 

p = ki + ka + . . . + k; (4.6) 

with some I > 2. The remaining m = n — I + 1 independent momenta will be denoted 
as qi, . . .qm- Analysis of the diagrams shows that m is equal to the number of the 
1-irreducible subgraphs of the form {v'vv . . that is, subgraphs with one field v' 

and arbitrary number of the fields v. This is easily checked for the examples given in 
figure 4: n = / = m = 1 for the diagram b and n = I = 2, m = 1 for the diagrams 
d and e|2 It is also clear that the momenta q^ here are exactly the same as in the 
original variables, while kj in (14.61) are certain linear combinations of p and the original 
momenta of the type k. 

Let us introduce the /-dimensional "vector" K = {ki, k2, . . . , h} (where /cj = |kj| > 
for alH = 1, . . . , /), the corresponding unit vector n = K./K, and pass to the spherical 
coordinates: 

/»oo /»oo /»oo /» 

/ dki... dki= dKK^-^ j dn. (4.7) 
Jo Jo Jo J 

The integrations over the directions of the vector n are restricted by the inequalities 

rij > 0. For d ^ oo, the scalar products can be dropped, so that 

p2 = (ki + k2 + . . . + kif c^kl + kl + ... + kf = K\ 

S This statement docs not apply to the graphs hke a and c in figure 4, which have no "harmful" 
{v'v')i-i-[ subgraph and involve no momenta of the type k; the trick is not needed for such diagrams. 
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The remaining variables can be made dimensionless by the replacement Qi —* Wi = Qi/ K 
{i = l,...,m = n — l + l). In the new set of variables n, Wi and K only the latter is 
dimensional. Thus the dependence of the integrand on K is found from the dimensional 
considerations as a known power, the integral over K factorizes out from the remaining 
integrals in the form 

POO 

4= / dKK^-'''^''+^^F{K^) (4.8) 
Jo 

and explicitly cancels the integral in the denominator of (14. 5p . This confirms the 
independence of the result on the choice of the weight function F{K'^). 

Thus we are left with some integrals over the dimensionless variables n and Wi of 
the general form 

/ du I dwi ... dwm V'(n, wi, . . . Wm) (4.9) 

Jni>0 J J 

with a known dimensionless function ip. They appear rather simple: all the two-loop 
integrals are calculated exactly (that is, without the expansion in the Laurent series 
in e) and are expressed in terms of F functions, trigonometric functions and simpler 
quantities. 

From the calculation procedure described above one can conclude that in the 
limit d oo, the momenta of all the propagators {vv)o (but not those of the mixed 
propagators {vv')o\) behave as if they were orthogonal to each other, including the 
case when one of such propagators corresponds to the external momentum p. For 
any given diagram it can be checked a posteriori that the condition that all such 
momenta be orthogonal defines the stationarity manifold in the integration region in 
the steep est- descent calculation for (i oo. However, it would be difficult to derive 
the orthogonality condition starting from the stationarity equations, especially in the 
original variables, where some of the propagators {vv)o correspond to linear combination 
of "pure" integration momenta. Thus the trick described above appears extremely 
useful. 

Let us illustrate the general scheme by the explicit calculation of the diagram d 
from figure 4. For this diagram n = I = 2 and m = 1 and, in agreement with the 
general statements made above, it involves two momenta of the type k in (14.61) and 
one momentum q. All these momenta are assigned to the lines {vv)q as shown in 
figure 5; the remaining lines correspond to their linear combinations. The dotted line 
denotes the weight function -F(j9^), where the former external momentum is represented 
as p = ki -|- k2. 

The integrand in the numerator of the expression (14. 5 p has the form 



(1 2 Ty^2\ 3 j 2 1 1 

-f^l^ ) ^ TZlTi \ U2(U2 , „2\ + (U2 , „2\(T^2 , .2^ + 



X 
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Figure 5. The new choice of the independent integration momenta in a two-loop 
diagram with the "harmful" subgraph. 



Here the first factor comes from tlie vertices (13.81) . tlie second one - from the time 
integrations and the third one - from the propagators {vv)q in (12.131) . After the angular 
integrations the integral in the numerator of the expression (14. 5 p takes on the form 

^3 roo POO POO 

J = - ,^o/o / dki / dk2 / dq / x 



128 (27r)2'^ io io ^ K^k.hqf^-^ 

+ ..2 , 2r.. , r...rl , (4-10) 



Passing to the new variables ki = Kui, q = Kw gives 

T - TT 

~ 128 (27r)2'^ °' 



where the integral 

POO 

/o= / dKK^-^'F{K^ 
Jo 



from (14. 8p cancels out the denominator in (14.51) and the remaining integral J has the 
form 

J = - I dn I dw i-^— Y < 2/ 2 I ^+ 



+ 7-^— ^T7^^ + r (4-11) 



'ni>0 JO 
1 

{ni + w^{l + w^) ' (1 + ^2) 

In the spherical coordinates ni = cos 6 and n2 = sin 6' the integral over the directions in 
(14. lip takes on the form 

Pn/2 

dn = dO. 

'ni>0 Jo 

The remaining integrations over w and 6 are performed explicitly and give the following 
final answer for the diagram in figure 5: 

J = r(.) r(-.) r(i - .) ( (ii^iWi^ + i2-e)ni-sn 

^' ^ ' ^ '\ 2r(2-3e) 4r(2-2£) J ^ ^ 

The explicit result (14.120 illustrates some general properties of the diagrams of the 
pair correlator (14.11) . In general, the diagrams contain poles at e = (manifestation of 



Renormalization group in the infinite- dimensional turbulence 



24 



the UV divergences) which are cancelled by the poles coming from the renormalization 
constants in the action fl2.17p in every order of the renormalized perturbation theory, 
so that the function R in (14.21) is finite at e ^ in every order of the expansion in u. 
In addition, the expression (14.121) involves the poles for some finite values of e, namely 
e = 1/2, 1, 3/2 and 2. All the other two-loop diagrams contain the poles at e = 1 and 
2 and some more diagrams contain the poles at £ = 1/2 and 3/2. Such singularities are 
manifestation of the IR divergences (that is, divergences coming from the region of small 
integration momenta) that occur in the diagrams if the calculation is performed in the 
"massless" (m = 0) model. When the IR cut-off m > in (12.91) is restored, the diagrams 
become convergent for all e > 0, and the IR problem is reformulated as the analysis 
of the singularities at m 0. The problem of such singularities in the stochastic NS 
equation with the stirring force ( 12. 9p was discussed in a number of studies; see e.g. 

m-m- 

The one-loop diagrams only contain the first-order poles at e = 1 and 2. As 
illustrated by the present two-loop calculation, in the higher-order diagrams the order 
of such poles increases and they are getting closer and closer to the origin e = 0. Thus the 
problem of the IR singularities lies beyond the scope of any given-order approximation 
of the RG, let alone the ordinary perturbation theory. For the correct analysis of the 
small-m behaviour, the standard RG techniques should be combined with additional 
methods - infrared perturbation theory or short- distance operator-product expansion 
[T0]-[T2]. For small values of e, this approach allows one to prove the existence of the 
limit m — >■ 0, which justifies the expansion in e of the scaling amplitude i?(l,M*) in 

The final result for the scaling function in (14. 2p has the form 

^^^'^)-2 + T6{2+H^-¥jhw + ---' ^'-''^ 

up to the terms of higher orders in u ^ e. Substituting the fixed-point coordinate (13.211) 
gives 

R{1, u^) = 1/2 + e/12 + (5/36 - TrVlOS)^^ + 0{e^). (4.14) 

This is the main result of the present section; it will be used below in the calculation 
of the Kolmogorov constant and the skewness factor. Although only the first two terms 
of the coordinate (I3.2ip are involved in the derivation of (I4.14p . the 0{e^) term will be 
needed to calculate the three terms of the whole expression (14. 2p . which contains the 

1 /3 

product uj -R(1,M*). In this sense, the two-loop calculation of the pair correlator is 
consistent with the three-loop calculation of the RG functions. 

II More physically, the IR singularities at e = 1/2 and 1 are related to the so-called sweeping effects 
and must disappear in the sum of all the diagrams for the equal-time correlator, which is a Galilean- 
invariant quantity; they do. The singularities at e = 3/2 and 2 do not disappear; they can be related 
to the fluctuations of the energy flux. These interesting issues will be discussed elsewhere. 
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5. Two-loop calculation of the Kolmogorov constant and skewness factor 

The Kolmogorov constant Ck can be defined as the dimensionless coefficient in the 
inertial-range asymptotic expression S'2(r) = Cxi^rY^^ for the second-order structure 
function, predicted by the classical phenomenological Kolmogorov-Obukhov theory 
[H EBjI^I Here £ is the average energy dissipation rate (cf. eq. (12.41) ) and the n-th 
order (longitudinal, equal-time) structure function is defined as 

Sn{r) = {[vrit,x + v) - Vrit,x)]''), Vr = {viri)/r, r = \r\. (5.1) 

Using this definition, the function 5*2 can be related to the momentum-space pair 
correlation function D{k) from ( 14. ip as follows: 

^ Dik) [1 - (kr)V(fcr)2] {1 - exp [i(kr)]} . (5.2) 

Alternatively, the Kolmogorov constant can be introduced through the 

phenomenological relation E(k) = C'j^£'^^^k~^^^, where the energy spectrum E{k) is 
related to the function (14.11) as E{k) = Sd{d — l)k'^^^D{k)/2. From the definitions one 
can derive the following relation between these two constants: 

3 • 2V^r(2/3)r(rf/2) ^, 

(rf + 2/3)1(^^/2 + 1/3)^^' ^^-^^ 

cf. [38] for d = 3. Using the exact relation S'3(r) = —12Sr/d{d + 2) that follows from 
the energy balance equation (see e.g. [H [3H] for d = 3), the constant Ck can be related 
to the inertial-range skewness factor S as follows: 

S = S^/ST = -\12 /d{d + 2)] C~k''\ (5.4) 

All these relations refer to real physical quantities in the inertial range, which in the 
stochastic model (12.11) . (12. 2p corresponds to e = 2 and m = in the random force 
correlator (12. 9p . 

Much work has been devoted to derivation of the Kolmogorov constant within the 
RG approach; see e.g. (TUl [HI EH SD 1121 113] and references therein. In order to 
obtain Ck, one usually combines the RG expression (14. 2 p for D{k) or the analogous 
expression for 5*2 by some relation between the physical parameter £ and the amplitude 
Dq in the random force correlator (12. 9p . In particular, in [H] the first-order term of 
the e expansion for the pair correlator was combined with the so-called eddy-damped 
quasinormal Markovian approximation for the energy transfer function, taken directly at 
the physical value e = 2. More elementary derivation, based on the exact relation (12.40 
between £ and the function df{k) from ( 12.90 was given in jllj; see also |5l|6]. In spite of 

% Now it is generally believed that the real exponent deviates slightly from 2/3 due to the phenomenon 
of intermittency; see the discussion in [1] and references therein. However, the situation is not absolutely 
clear: some researchers argued that the existing experimental data are consistent with the "2/3 law" 
and the observed disagreement is due to the corrections (and not deviations) vanishing in the limit 
of infinite Reynolds number; see e.g. [40j . Now we are interested in the large-d limit, in which 
intermittency is expected to get weaker or to disappear; so we will accept the "2/3 law," which is 
also internally consistent within the calculation based on the e expansion. 
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a reasonable agreement with the experiment, such derivations are not satisfactory from 
the theoretical viewpoints. Their common drawback is that any relation between S and 
Do is unambiguous only in the limit e ^ 2 (see equation f l2.8p ). so that the coefficients 
of the e expansions for Ck can be made arbitrary; see e.g. the discussion in [32] and 
Sec. 2.10 of [6]. 

This ambiguity is a consequence of the fact that the notion itself of the Kolmogorov 
constant has no unique extension to the nonphysical range < e < 2. Furthermore, 
the experience from the RG theory of critical behaviour suggests that well-defined 
e expansions can be written for universal quantities, such as critical exponents or 
normalized scaling functions, which do not involve bare parameters. The constant Ck 
extended to the range < e < 2 as in Refs. [TT|, HI], |12] involves a bare parameter, Dq, 
and hence is not universal in this sense. 

To circumvent these difficulties, an alternative derivation was proposed in that 
does not involve any relation between Dq and S, relates Ck to a universal (in the sense of 
the theory of critical behaviour) quantity, and leads (for c? = 3) to reasonable agreement 
with experimental data. Below we adopt this derivation to our case, d ^ oo. 

Consider the ratio 

Q{e) = VrS2{r)/\Ss{r)\'/' = V^S^ir) / {-S^{r)f'' (5.5) 

with T>r = rd/dr. We shall see below that it is universal and can be calculated in the 
form of a well-defined e expansion. On the other hand, its value at e = 2 determines 
the Kolmogorov constant and the skewness factor through the exact relations 

Ck = [3g(2)/2] [U/d{d + 2)f^ , S = - [3g(2)/2]-'/' (5.6) 

which follow from the definitions, relation (15. 4p and the identity Vrr^ = 6r^ for any 6. 

From the RG equations one derives the analogs of the representation (14.21) for the 
structure functions entering the ratio (15. 5p : 

Ss{r) = Dor-'^'-fsie), V^S^ir) = dI'\-'''^ h{e). (5.7) 

Here m = 0, /xr > 1, < e < 2 and = 1 - 2£:/3, cf. (g^]). The operation 
introduced in (15. 5p "kills" the constant contribution (f ^) in 5*2 which does not exist 
without an UV cutoff for e < 3/2 (see below); in 5*3 such contribution is absent. 

It follows from (15. 7p that the amplitude Dq disappears from the ratio Q{s) = 
/2/(-/3)^/^ and the latter can be calculated in the form Q{e) = e^^^p{e), where p{e) 
is a power series in e. Our results for the fixed point from section [3] and for the pair 
correlator from section H] allow us to find the first three terms of p{e). In this sense, one 
can speak about the third-order approximation for the Kolmogorov constant (previous 
attempts have mostly been confined with the first order, with the second-order exception 
of Ref. |10j). To avoid possible misunderstanding, we stress again that we did not 
intend to extend the definition of the physical quantities Ck and S to the whole interval 
< e < 2 and to construct their e expansions from the known expansion for Q{e). 
Instead, the latter is used to give the value of Q(2), which, in its turn, determines Ck 
and S through the relations (15. 6p that make sense only for the real value e = 2. 
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Applying the operation to the expression fl5.2p gives: 

VrS2{r) = 2 j D{k) [I - (kr)V(fcr)2] (kr) sin(kr). (5.8) 

In order to obtain the inertial-range form of the function VrS2{r) it is sufficient to 
substitute the asymptotic expression (14. 2 p into (15. 8p . A straightforward calculation 
gives: 

- w4'/2;2y3) «^ ^«"(^/^)-^" (^■^) 

with the amplitude /2(1,m*) from f l4.2p . It is important here that the integral in (15. Sp 
with the function (14. 2p exists for all < e < 2. It is the operation "D^ that ensures its 
convergence: the original integral (15. 2p would be UV divergent for < e < 3/2. 

The needed terms of the e expansion for /s can be obtained from the direct 
perturbative calculation, similar to the calculation of the pair correlator in section HJ 
but it is more convenient to use the exact expression 

which follows from the energy balance equation and generalizes the well-known 
Kolmogorov's "4/3 law" [Il[38] to general d and e. Thus from Eqs. (15. 9p and (I5.10p for 
the ratio Q{e) one obtains 

Q{e) = [4{d - 1)m,/9]'/^ A{e) R{1, u,), (5.11) 

where we excluded in favor of = g^Sd] the coefficient A{e) is given by 

._ r(2-2e/3)rV3(rf/2)r^/3(rf/2 + ,) 

r(t//2 + 2e/3)P/3(2-5) • ^^■'^> 
Using the Stirling formula for the F functions one can show that 

r(<i/2 + 2e/3) 
so that in the limit d ^ oo one obtains 

^(-) = %]2-e) ='^U'- "'/^^ + ^'-''^ 
Substituting the expressions (I3.2ip . (I4.14p and (I5.13P into (15.1 ip gives 

(4ec/)i/3 e ^ f 49 



«'^) = ^ + 18 ni62 - 27 j ^ + )| ■ f^-"' 
Now the Kolmogorov constant Ck and the skewness factor S are obtained from the 
expressions (15. 6p . with the replacement d{d + 2) d'^ in the first one. Thus in the 
leading order of the large-rf asymptotic behaviour one obtains Ck ex. 1/d and S oc 1/ d^^"^. 
From the relation (15. 3p it then follows C'j^ oc d^^^, in agreement with the earlier results 
derived within the direct interaction approximation [2^ and within the RG approach 

EE]. 
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Although these results refer to the large-c? limit, one can try to use them as some 
approximation to the real three-dimensional case. Substituting the third-order result 
f l5.14p into the second relation from (15.61) and setting e = 2 and d = 3 gives S ~ —0.73, 
while the experimental value recommended in [38j is 5 ~ —0.28. Of course, one should 
not have expected a better agreement. Surprisingly enough, the situation appears 
much better for the Kolmogorov constant. Substituting the first, second and third 
approximations from fl5.14p into the first relation from (15. 6p . one obtains for 5 = 2 and 
d = 3 the following results: 

dj^^ ^ 1.75, CP ^ 1.94, cP ^ 1.50, 

all of them in a reasonable agreement with the experimental estimate Ck ~ 1-9 
recommended in [3H]. One can speculate that the skewness factor which is related 
to the odd-order function 5*3 and is a measure of anisotropy in the distribution of the 
velocity field is more sensitive to the spatial dimension than the Kolmogorov constant, 
related to the even-order function 82- 



6. Conclusion 



We have accomplished the complete three-loop calculation of the main ingredients of 
the RG analysis of the d-dimensional stochastic Navier-Stokes equation (12. ip . (12.20 in 
the limit d 00: the renormalization constant Z,y, the (3 function, the coordinate of 
the fixed point and the ultraviolet (UV) correction exponent u; these results are 
summarized in equations (I3.19p - (l3.22p . We have also calculated in the large-c? limit 
the correlation function of the velocity in the third order of the e expansion (two-loop 
approximation) . 

These results allowed us to derive the third-order answers for two interesting 
physical quantities: the Kolmogorov constant Ck in the spectrum of turbulent energy 
and the inertial-range skewness factor S. The former is in a reasonable agreement with 
the existing experimental data for the real three-dimensional turbulence. 

The successful analytic calculation in the third order appeared feasible due to 
drastic simplifications that we found in the large-c? limit: in particular, most of the 
Feynman diagrams for the Green function in that limit vanish, and the remaining ones 
are reduced to relatively simple (and analytically calculable) integrals. Some more 
efforts and tricks were required to derive the third-order results for the pair correlation 
function of the velocity. 

Although we did not succeed in finding the exact solution for c? = 00, the 
simplifications that occur in the calculations in that limit and the simple form of the 
obtained results suggest that this is not an impossible task. Based on the three-loop 
answers, we proposed some hypothetical expressions beyond the e expansion which show 
what, in principle, such exact results may look like. 

We believe that the present results, as well as the calculational techniques developed 
in their derivation, will be useful in the further attempts of the construction of the 
systematical 1/d expansion for the fully developed turbulence. 
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